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ABSTRACT 

We present results from a Chandra observation of the NGC 346 cluster. This cluster contains 
numerous massive stars and is responsible for the ionization of N66, the most luminous Hn 
region and the largest star formation region in the SMC. In this first paper, we will focus on the 
characteristics of the main objects of the field. The NGC 346 cluster itself shows only relatively 
faint X-ray emission (with L^ abs ~ 1.5 x 10 34 erg s _1 ), tightly correlated with the core of the 
cluster. In the field also lies HD 5980, a LBV star in a binary (or possibly a triple system) that 
is detected for the first time at X-ray energies. The star is X-ray bright, with an unabsorbed 
luminosity of l™ ,ahs ~ 1.7 x 10 34 erg s _1 , but needs to be monitored further to investigate its 
X-ray variability over a complete 19 d orbital cycle. The high X-ray luminosity may be associated 
either with colliding winds in the binary system or with the 1994 eruption. HD 5980 is surrounded 
by a region of diffuse X-ray emission, which is a supernova remnant. While it may be only a 
chance alignment with HD5980, such a spatial coincidence may indicate that the remnant is 
indeed related to this peculiar massive star. 

Subject headings: (galaxies:) Magellanic Clouds-X-rays: individual (NGC 346, HD 5980) 
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1. Introduction 

The Small Magellanic Cloud (SMC) is an irreg- 
ular galaxy at a distance of 59 kpc (Mathewson, 
Ford, & Visvanathan 1986) that forms a pair 
with the Large Magellanic Cloud (LMC). Both 
are satellites of our own Galaxy. The interstellar 
extinction towards the Magellanic Clouds is low, 
allowing studies of the X-ray sources to be under- 
taken with clarity. 

Previous X-ray observations of the SMC have 
been made with the Einstein Observatory, ASCA 
and ROSAT. These observations included surveys 
of the point source population (Kahabka et al. 
1999; Schmidtke et al. 1999; Haberl et al. 2000; 
Sasaki, Haberl, & Pietsch 2000; Yokogawa et al. 
2000) and in particular studies of the properties 
of the X-ray binary population (Haberl & Sasaki 
2000). The recent launch of Chandra, however, 
provides an opportunity to study the SMC with a 
far greater sensitivity and spatial resolution than 
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ever before. 

We have obtained a Chandra observation of 
the giant Hn region N66 (Henize 1956), the 
largest star formation region in the SMC, in 
order to study the young cluster NGC346 and 
its interaction with the surrounding interstellar 
medium. This cluster contains numerous massive 
stars (Massey, Parker & Garmany 1989), of spec- 
tral types as early as 02 (Walborn et al. 2002). 
The large number of massive stars is not the only 
feature of interest in this field. On the outskirts of 
NGC 346 lies the remarkable star HD 5980, which 
underwent a Luminous Blue Variable (LBV)-type 
eruption in 1994. This massive binary (or triple?) 
system has been monitored and analysed for its 
varying spectral and photometric properties in 
visible and UV wavelengths since the early 80's 
(Koenigsberger et al. 2000; Sterken & Breysacher 
1997, and references therein). Previous X-ray ob- 
servations of NGC 346 (Inoue, Koyama, & Tanaka 
1983; Habcrl et al. 2000) detected a bright ex- 
tended source around HD 5980 which has been 
attributed to a supernova remnant (SNR); how- 
ever, the crude instrumental resolution prohibited 
unambiguous detections of a point source associ- 
ated with HD5980. In addition to HD5980, the 
young massive stars of NGC 346 and their inter- 
acting winds are likely to produce X-ray emission, 
but have not been detected previously. The sen- 
sitive, high-resolution Chandra observation of the 
NGC 346 field thus provides an excellent oppor- 
tunity to study a variety of phenomena involving 
emission at X-ray energies. 

In this paper, we will describe in § 2 the ob- 
servations used in this study, then discuss the 
available data on NGC 346, HD 5980, and the ex- 
tended emission in § 3, 4, and 5, respectively. 
Finally, conclusions are given in § 6. The second 
part of this work (Naze et al., paper II) will de- 
scribe the properties of the other sources detected 
in the field. 



2. The Observations and Data Analysis 
2.1. X-ray Observations 

NGC 346 was observed with Chandra for the 
XMEGA 9 consortium on 2001 May 15-16 for 
100 ks (98.671 ks effective, ObsID = 1881, 
JD~2 452 045.2d). The data were recorded by 
the ACIS-I0, II, 12, 13, S2 and S3 CCD chips 
maintained at a temperature of — 120°C. The data 
were taken with a frame time of 3.241s in the faint 
mode. The exposure was centered on the cluster, 
with HD5980 lying 1.9' to the north-cast of the 
aimpoint on ACIS-I3. Our faintest sources have 
fluxes of about 2 x 10 32 erg s _1 (see paper II), 
assuming a distance of 59 kpc for the SMC. 

We excluded bad pixels from the analysis, us- 
ing the customary bad-pixel file provided by the 
Chandra X-ray Center (CXC) for this particular 
observation. We have searched the data for back- 
ground flares, which are known to affect Chandra 
data, by examining the lightcurve of the total 
count rate, but no flares were found. Event Pulse 
Invariant (PI) values and photon energies were 
determined using the FITS Embedded File (FEF) 
acisD2000-01-29fef_piN0001.fits. 

For long exposures, removing the afterglow 
events can adversely affect the science analysis 
(underestimation of the fluxes, alteration of the 
spectra and so on) 10 . We thus computed a new 
level 2 events file by filtering the level 1 file and 
keeping the events with ASCA grades of 0, 2, 3, 
4 and 6, but without applying a status=0 filter. 
Throughout this paper, we will use this new file 
for all scientific analysis. 

Further analysis was performed using the CIAO 
v2.1.2 software provided by the CXC and also 
with the FTOOLS tasks. The spectra were anal- 
ysed and fitted within XSPEC vl 1.0.1 (Arnaud 
1996). In Fig. 1, we show a 3 color image of the 
Chandra data of the NGC 346 cluster, HD 5980 
and its surroundings. We will discuss the features 
of this image further in later sections. 



9 http://lheawww. gsfc.nasa.gov/users/corcoran/xmega/xmega.html 
°see caveat on http://asc.harvard.edu/ciao/caveats/acis_cray.html 
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We have investigated problems due to Charge 
Transfer Inefficiency (CTI). To remove these CTI 
effects, we used the algorithm of Townsley et al. 
(2000). Checking several sources on different chips 
and/or at different positions on each CCD, we con- 
cluded that the impact of the CTI was small and 
that it did not change significantly the spectral 
parameters found by using non-corrected data. 
When a difference was apparently found, we dis- 
covered that the best fits found on corrected and 
non-corrected data gave similar \ 2 ■ As the work 
of Townsley et al. (2000)is not yet an official 
CXC product 11 , we chose to present here only the 
results of the non-corrected data. 

2.2. HST WFPC2 Images 

Two HST WFPC2 observations of the N66 
/ NGC346 region in Ha are available from the 
STScI archives. They were taken on 1998 Sept. 27 
and 2000 August 7 for the programs 6540 and 
8196, respectively. HD 5980 is centered on the PC 
in the observations from program 6540 (see Fig. 
2), while program 8196 shows an area situated fur- 
ther south. The observations were made through 
the F656iV filter for 2070s (Prog. # 6540) and 
2400s (Prog. # 8196). This filter, centered at 
6563.7 A with a FWHM of 21.4A, includes the 
Ha line and the neighboring [N n] lines. 

The calibrated WFPC2 images were produced 
by the standard HST pipeline. We processed them 
further with IRAF and STSDAS routines. The 
images of each program were combined to remove 
cosmic rays and to produce a total-exposure map. 
No obvious correlation between X-rays features 
and the Ha emission is seen in the HST data. 
Fig. 2 will be discussed more extensively in § 5. 

3. The NGC 346 Star Cluster 

The NGC 346 cluster has a large number of 
massive stars, containing almost 50% of all early- 
type O stars in the SMC (Massey, Parker & Gar- 
many 1989). It is responsible for the ionization 
of N66, the most luminous Hn region in the SMC 
(Ye, Turtle, & Kennicutt 1991). However, the 

11 http://www.astro.psu.edu/users/townsley/cti/ 



NGC 346 cluster has not been detected in X- 
rays until this Chandra observation. Even though 
NGC 346 lies near the joint of the four ACIS-I 
CCD chips and part of the cluster is lost in the 
gaps between CCDs, the X-ray emission from the 
core of NGC 346 is unambiguously detected in the 
corner of the ACIS-I3 chip. The X-ray emission 
appears extended with multiple peaks correspond- 
ing to the bright blue stars MPG 396, 417, 435, 
470, and 476 (Fig. 3). As only ~300 counts are de- 
tected and the X-ray peaks are not well-resolved, 
we will analyze only the overall emission from the 
cluster. 

The total count rate of this emission in the 
0.3 - 10.0 keV band is 3.04 ± 0.36 x 10~ 3 cnts s" 1 . 
The spectrum of the NGC 346 cluster was ex- 
tracted from a square aperture of ~45" x 45" with 
a carefully chosen background region, as close as 
possible to the cluster. As we regard the clus- 
ter X-ray emission to be extended, we have used 
the calcrmf and calcarf tools to generate the ap- 
propriate response files. The best fit (x 2 =0.90, 
Ndof=50, Z = O.IZq) to this spectrum was an 
absorbed mekal model (Kaastra 1992) with the 
following properties (see Fig. 4): an absorption 
column density N(H) of 0.42g;ff x 10 22 cm~ 2 and 
a temperature kT of 1.01 2 ,; 2 ; 3 keV. The luminosity 
in the 0.3 - 10.0 keV band was 5.5 x 10 33 erg s" 1 , 
which corresponds to an unabsorbed luminosity 
of 1.5 x 10 34 erg s _1 . The absorption column is 
comparable to the value expected for NGC 346 
(see paper II) , but the temperature is one of the 
lowest we found, except for the extended source 
around HD 5980. 

In the region of X-ray emission in NGC 346, 
30 blue stars were detected by Massey, Parker & 
Garmany (1989): MPG 435, 342, 470, 368, 476, 
396, 487, 417, 370, 471, 467, 495, 451, 330, 455, 
454, 500, 445, 468, 375, 508, 395, 439, 371, 485, 
561, 374, 366, 557, 400 (by increasing magnitude). 
Of these, 16 have known spectral types that we 
can use to convert magnitudes to bolometric lumi- 
nosities. The bolometric correction factors were 
taken from Humphreys & McElroy (1984). We 
can then compute the expected X-ray luminosi- 
ties using log{L u ^ abs ) = 1.131og(L B oz,) - 11-89 
(Berghofer et al. 1997). If the X-rays were com- 
ing only from these 16 stars, we should expect a 
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total unabsorbed luminosity of 2 x 10 33 erg s _1 . 
This represents only 13% of the detected luminos- 
ity. Several hypotheses could explain this discrep- 
ancy: a metallicity effect (Bcrghofcr's relation was 
determined for Galactic stars but SMC stars have 
weaker winds) ; X-ray emission from the other (un- 
classified) stars in this region, especially from the 
low- mass stars (see e.g. Getman ct al. 2002); ad- 
ditional emission linked to interactions in binary 
systems (such as colliding winds); or an extended 
region of hot gas as the winds from the individual 
stars combine to form a cluster wind (Ozernoy, 
Genzel, & Usov 1997; Canto, Raga, & Rodriguez 
2000). 



Chandra has also observed X-ray emission from 
other very young (age < 5 Myr) stellar clusters, 
comparable to NGC 346. The Arches cluster in the 
Galactic Center region has been observed to have 
emission from several different regions, with a to- 
tal X-ray luminosity of £™ abs ~ 5 x 10 35 erg s _1 
(Yusef-Zadeh et al. 2002). NGC 3603 has also 
been observed with Chandra (Moffat et al. 2002) 
and the total unabsorbed cluster luminosity is 
L unabs _ x x 1Q 35 erg g -i (though in this case it 
is clear that a substantial fraction of the emission 
is from point sources in the cluster). For both 
the Arches and NGC 3603 clusters the fitted tem- 
perature of the X-ray emission is higher than the 
kT ~ 1 kcV value for NGC 346. 

It is important to note that the calculations 
presented by Canto, Raga, & Rodriguez (2000) 
and Raga et al. (2001) for the Arches cluster as- 
sume that there are about 60 stellar sources hav- 
ing mass-loss rates of 10 -4 M© yr _1 . If we take 
the list of the bluest stars in NGC 346 given by 
Massey, Parker & Garmany (1989), the first 60 
sources have spectral types in the range O3-B0. 
Assuming mass-loss rates for 8 SMC stars given 
by Puis et al. (1996) to be typical, then all stars 
in NGC 346 (excluding HD5980) have mass-loss 
rates smaller than 10~ 4 M Q yr _1 . Only MPG 
435 (04III(n)(f)) has a mass loss rate as high as 
9xl0 -5 Mq yr _1 . However, the majority of these 
60 stars are late O-type Main Sequence stars or 
cooler, which must have much smaller mass-loss 
rates (10 -6 M© yr _1 or less), and they are more 
widely distributed than in Arches. Thus, it is not 



surprising that the X-ray luminosity is lower than 
in Arches, and, in fact, it may be too large to be 
explained solely with the cluster wind scenario. 

4. The LBV HD 5980 

Among all stars in NGC 346, HD 5980 is unique 
in its spectral variations in the last two decades, 
and warrants further examination. This eclips- 
ing binary (or triple system?) was classified 
as WN+OB before 1980, but its spectral type 
changed to WN3+WN4 in 1980-1983, and to 
WN6, with no trace of the companion in 1992. 
In 1994, this star underwent a LBV-type erup- 
tion, and presented at the same time a WN11 
type. The eruption has now settled down and the 
spectrum is returning to its pre-eruption state (for 
more details see Koenigsberger et al. 2000, and 
references therein). We will label the 1994 emptor 
as 'star A' and its companion as 'star B\ 

Chandra is the first X-ray telescope to detect 
this peculiar system, since the low spatial resolu- 
tion of previous X-ray observatories did not allow 
the distinction of HD 5980 from the extended emis- 
sion surrounding it. However, even with a 100 ks 
exposure, the data possess a rather low signal-to- 
noise ratio. The Chandra ACIS-I count rate of 
HD 5980 is only 2.98 ± 0.19 x 10~ 3 cnts s" 1 in the 
0.3 - 10.0 keV band or about 300 counts in the 
total observation. This limits our ability to sensi- 
tively determine the spectral shape of the emission 
and to look for source variability. On the other 
hand, the relative faintness of the source means 
that photon pileup will not affect our analysis to 
any great degree. 

We have extracted a spectrum of HD 5980 using 
the CIAO tool psextract. An annular background 
region around the source (but sufficiently far from 
it) was chosen, in order to eliminate contamina- 
tion from the surrounding extended emission (see 
§ 5). This spectrum is shown in Fig. 5a. The 
best fit model to the spectrum has the follow- 
ing parameters: N(H) = 0.22^; 3 | x 10 22 cm" 2 
and kT = 7.04| 3 3 f 5 keV for an absorbed mekal 
model and N(H) = 0.28g;f| x 10 22 cm" 2 and 
r=1.74};53 for an absorbed power-law. These ab- 
sorption columns are consistent with the value 
expected for NGC 346 (see paper II) but part of 
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these columns are probably due to wind absorp- 
tion and/or absorption from the 1994 ejecta. The 
derived observed X-ray luminosity of HD 5980 is 
L x = 1.3 x 10 34 erg s" 1 in the 0.3 - 10.0 kcV band 
for both models. Using the normalisation factor 
of the mekal model, we found a volume emission 
measure of ~10 57 cm -3 for this source. 

We can compare HD 5980 with other WR stars 
detected in X-rays. The unabsorbed X-ray lu- 
minosity of HD5980 is 1.7 x 10 34 erg s" 1 in the 
0.3 - 10.0 keV energy range and 9 x 10 33 erg s" 1 
in the ROSAT range, i.e. 0.2 - 2.4 kcV. This 
places HD 5980 amongst the X-ray brightest single 
WN stars (Wessolowski 1996) and the brightest 
WR+OB binaries (Pollock 2002) of the Galaxy. 

Two factors could explain this high X-ray lu- 
minosity. First, the fast wind from the post- 
cruptivc phases (from 1300 kms -1 in 1994 to 
-2000 kms" 1 in 2000) should now collide with 
the slow wind (ejected with a velocity as low as 
200 kins -1 during the eruption, see e.g. Koenigs- 
berger et al. 2000), which will increase the post- 
eruptive X-ray luminosity. X-ray emission from 
colliding ejecta has already been observed for an- 
other LBV, i.e. J] Carinae. An clliptically shaped 
extended emission of 40"x70"(i.e. 0.4 pcx0.8 pc) 
in size surrounds the star in the X-ray domain 
(Seward et al. 2001). It is correlated to the high 
velocity ejecta of the Homunculus nebula (Weis, 
Duschl, & Bomans 2001), which was created af- 
ter the last great eruption. For HD 5980, however, 
the ejecta from the 1994 eruption has not had 
sufficient time to form a detectable LBV nebula 
around the star, and no 'Homunculus-like' nebula 
from a previous eruption is visible around HD 5980 
(see Fig. 2): any X-ray emission from the colliding 
ejecta should thus still be blended with the stellar 
emission from HD 5980 and that may contribute to 
explain the high X-ray luminosity. Unfortunately, 
no X-ray detection of HD5980 before the 1994 
eruption is available. A ROSAT all-sky survey im- 
age taken in 1990 (exposure number 933001) pro- 
vides only an upper limit of Lx — 10 36 erg s _1 at 
the position of HD 5980. Wang & Wu (1992) gave 
a 3ct upper limit of 2xl0 34 erg s" 1 for all SMC 
WR stars, except AB7. These limits are compat- 
ible with the detected luminosity of HD 5980 and 
do not allow us to quantitatively compare the prc- 



and post-eruption luminosities to confirm the pre- 
dicted luminosity enhancement. 

The fact that HD5980 is a close binary sys- 
tem containing two very massive stars, suggests 
that colliding winds may provide another source 
for the observed X-rays. We can thus estimate 
the expected X-ray luminosity. For the sys- 
tem, we assume that the stellar winds of the 
stars have returned to their pre-eruption values 
and we assume that M,4=1.4x 10~ 5 M yr -1 , 
Woo (yl)=2500 km s" 1 for the O-star, and M B = 
2xl0~ 5 M yr' 1 , v oc (B)=1700 km s" 1 for the 
WN Wolf-Rayet star (Moffat et al. 1998). Thus 
the momenta of each star's wind are very similar, 
and the total wind kinetic energy is 5 x 10 37 erg 
s _1 . These wind parameters lead to a momen- 
tum ratio of r/ = (M B w 00 (B))/(M j4 i; 00 (A)) = 0.97 
(Usov 1992) and the wind-wind collision shock 
will lie halfway between the stars. In this config- 
uration we expect about 1/6 of the wind kinetic 
energy to pass perpendicularly through the two 
shocks and be thermaliscd (Stevens, Blondin, & 
Pollock 1992), and because the systems are rela- 
tively close much of this energy will be radiated 12 . 

Consequently, we would expect the X-ray lumi- 
nosity of HD 5980 to be — 10 36 erg s" 1 or more, 
rather than the observed value of ~ 10 34 erg s _1 . 
The discrepancy between the predicted and ob- 
served X-ray luminosities could be due to a range 
of factors, e.g. the winds may not be at the pre- 
eruption levels - lower mass-loss rates or lower 
wind velocities will translate to a lower X-ray lu- 
minosity. Radiative braking, wind clumping, or 
the eclipse of the colliding wind region by the 
thick wind of star B could also reduce the lumi- 
nosity. Alternatively, the wind momenta may not 
be so nearly equal, as a low value of r\ also tends 
to lower the luminosity. 

Using the ephemeris of Sterken & Breysacher 
(1997), we have computed a phase of (f> = 0.24 — 
0.30 for our Chandra observation: this is close to 
the eclipse of star A by star B (</> ec ;=0.36). If the 
eccentricity is e — 0.3, a simple adiabatic model 



In terms of the cooling parameter \ denned by Stevens, 
Blondin, & Pollock (1992), for the WN wind % < 1 and 
for the O-star wind \ ~ 1. 
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predicts that the X-ray luminosity will vary in- 
versely with separation. We may expect a change 
in the intrinsic X-ray luminosity of a factor of ~ 2 
through the orbit and by ~ 10% during the Chan- 
dra observation. We have detected a variation 
of apparently larger amplitude (factor ~2) in the 
count rate of HD 5980 (see Fig. 5b). This still sug- 
gests we are seeing orbital variability associated 
with colliding winds, since the additional effects 
of changing absorption 13 on the observed luminos- 
ity can either reduce or enhance this variability. 
This observed variability is much more likely to 
be due to colliding winds than wind-blown bubble 
type emission, where no short timescale variability 
would be expected. The high fitted temperature 
of HD 5980 kT ~ 7 keV also suggests that we are 
likely seeing colliding wind emission, as this tem- 
perature corresponds to a shock velocity of <~ 2500 
km s _1 , comparable to the wind speed of star A 
(voc(A) above). 

The detailed characteristics of the X-ray prop- 
erties of HD 5980 need to be studied further, in 
conjunction with detailed colliding wind models 
(c.f. the case of r\ Carinae; Pittard & Corcoran 
2002). An XMM-Newton observation is scheduled 
at phases <fi — 0.09 — 0.10, just after periastron: 
the comparison between these two datasets may 
enable us to detect further the signature of the col- 
liding wind region, and constrain it more precisely. 

Finally, we note that the unabsorbed X-ray 
luminosity of the central source in 77 Carinae, 
if we assume a distance of 2.3 kpc, varies be- 
tween (0.6 and 2.5) xlO 35 erg s _1 (Ishibashi et al. 
1999). Such luminosities are much higher than for 
HD5980. 

5. The Extended Emission around HD 5980 

Around the position of HD 5980 lies a region 
of bright extended X-ray emission. It was first 
detected by the Einstein Observatory (source IKT 
18 in Inoue, Koyama, & Tanaka 1983), and sub- 
sequently observed by ROSAT (source [HFP2000] 

The colliding wind region is viewed alternatively through 
the atmosphere of stars A and B. The absorption thus 
changes with the orbital phase. 



148 14 in Haberl et al. 2000). The Chandra image 
reveals the extended emission in much more detail 
than previous X-ray observatories. The overall 
shape of the emission region is more or less rect- 
angular, with an extension to the northeast. It 
contains a few bright or dark arcs, but apart from 
these, the brightness is rather uniform, with no ob- 
vious limb-brightening. It has a size of 130" x 100", 
i.e. 37pcx29pc at the SMC's distance. HD 5980 
lies towards the top center of the emission region. 
An X-ray bright filament extends from 9" west to 
23" south of HD 5980. An X-ray dark feature ap- 
pears some 30" at the southeast of HD 5980. The 
total count rate of this source in the 0.3 — 10.0 keV 
energy range is 7.25 ± 0.17 x 10~ 2 cnts s _1 . It is 
the softest source present in the field. 

The spectrum of this extended X-ray emis- 
sion was extracted in a rectangular aperture of 
~160"xll0" with a carefully chosen background 
region, located as close as possible to the source 
and in a region on the ACIS-I1 CCD where there 
are no point sources. A circular region of di- 
ameter ~5" containing HD 5980 was removed 
from the extraction. As the X-ray emission is 
extended the reponse matrices were calculated us- 
ing the calcrmf and calcarf tools. The best fit 
(X^=l-03, A^ /=244) to this spectrum was an 
absorbed mekal model with the following proper- 
ties (see Fig. 6): N(H) = 0.12g ^ x 10 22 cuT 2 , 
kT = 0.66g;g| keV and Z = 0.17%jlZ Q . The 
observed flux in the 0.3 — 10.0 keV band was 
3.35 x 10 -13 erg cm ~ 2 s _1 , i.e. an observed lu- 
minosity of 1.40 x 10 35 erg s _1 . The low value 
of the absorption column, compared to HD 5980 
and the cluster, suggests that this X-ray source 
lies between NGC 346 and the observer, but still 
in the SMC (see paper II). 

Using the normalisation factor of the mekal 
model, we found a volume emission measure 
/ n e nndV of ~3xl0 58 cm~ 3 . Considering a con- 
stant density for the hot gas, a sperical geometry 
of diameter ~33 pc for the extended emission and 
assuming a pure H composition, the density of the 
hot gas is roughly 0.2 cm~ 3 and the total mass of 
the hot gas is ^100 M Q . 

This ROSAT source may encompass some of the X-ray 
emission from the cluster. 
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Using the available data, a deeper analysis was 
conducted, searching for the existence of a temper- 
ature gradient throughout the source area, with 
color and temperature maps. For color maps, 
we generated a set of narrow-band images, e.g. 
0.3 - 0.9 keV, 0.9 - 1.2 kcV and 1.2 - 2.0 keV, 
that enable us to distinguish the harder compo- 
nents from the softer parts. In addition, temper- 
ature maps were constructed from spectral fits 
in small regions of the SNR. The only obvious 
trend in both the color and temperature maps is 
the softness of the northeastern extension of the 
SNR. However, it is not possible to draw any firm 
conclusions regarding the presence of any temper- 
ature gradient. 

We have also constructed images of the ex- 
tended emission in narrow energy bands, each 
correspond to a specific ion. The energy bands 
used for each ion are shown in Fig. 7. The data 
in each band were binned to obtain 4.9" x 4.9" 
pixels (see Fig. 7). In contrast to N132D (Bchar 
et al. 2001), no significative differences between 
the morphology of highly ionized species (Mg 10+ , 
Ne 9+ , Si 12+ , Fe 19+ ) and lower ionization species 
were detected. 

5.1. The Nature of the Extended Emission 

Since its discovery the extended X-ray emis- 
sion has been attributed to a supernova remnant 
(SNR). Further evidence was subsequently ob- 
tained to support this hypothesis; for example, 
Ye, Turtle, & Kennicutt (1991) found a non- 
thermal radio source at this position that they 
called SNR 0057 - 7226. Using the lowest con- 
tours of Ye, Turtle, & Kennicutt (1991), the 
radio source is 56 pcx52 pc, slightly larger than 
the X-ray source. 

Moreover, evidence of high velocity motions in 
this region were observed in the visible and UV 
ranges. Using echelle spectra centered on Ha, 
Chu & Kennicutt (1988) found clumps moving 
at vlsr ~ 300 km s -1 , rcdshifted from the main 
quiescent component by + 170kms _1 : only the 
receding side of the expanding object is seen, sug- 
gesting a low density where the approaching side 



arises. UV analyses with IUE (de Boer & Savage 
1980), HSTSTIS (Koenigsberger et al. 2001) and 
FUSE (Hoopcs et al. 2001) have also confirmed 
the presence of an expanding structure. While de 
Boer & Savage (1980) and Hoopes et al. (2001) 
found only absorptions at Vi5fl=+300kms _1 , 
Koenigsberger et al. (2001) detected both ex- 
panding sides of the object, with components at 
vlsr = +21, +52, +300, +331 kms -1 . However, 
the absorptions at +21 and +52 kms -1 were 
small, and need a future confirmation. The pres- 
ence of components from both approaching and 
receding sides of the expanding structure has now 
been confirmed by Danforth et al. (2002) from 
deep echelle spectroscopic integrations at optical 
wavelengths of the N66 region. Absorptions at 
any of these velocities are not seen in the spectra 
of Sk 80, a close neighbor of HD 5980 situated on 
the edge of the X-ray source. A fast expanding 
structure close to the position of HD 5980 is thus 
present. 

Unfortunately, even the high-resolution HST 
WFPC2 images of the close neighborhood of 
HD 5980 do not show any clear nebula associated 
with the X-ray extended emission (see Fig. 2). A 
lot of filamentary structures - generally indicative 
of a SNR (Chen et al. 2000) - are seen throughout 
the field, but are not limited to the exact location 
of the X-ray source. There is thus no obvious Ha 
feature directly correlated with the extended X- 
ray emission. 

Considering all the evidence (non-thermal emis- 
sion, high velocity expanding shell etc) , the X-ray 
/ radio source should thus be regarded as a SNR 
located in front of HD 5980 but still belonging 
to the SMC. The SNR hypothesis is further sup- 
ported by the comparison of the size of this fea- 
ture to its X-ray and radio fluxes (Mathewson et 
al. 1983). 

As HD 5980 is projected more or less at the cen- 
ter of the diffuse X-ray emission and it underwent 
a LBV-type eruption in 1994 (which is unlikely to 
be the only one undergone by the system), it is 
tempting to associate the diffuse X-ray emission 
with HD 5980. To investigate whether this asso- 
ciation is likely, we first compare HD 5980 with 
X) Carinae, the most well-known LBV that have 
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gone through multiple eruptions. 77 Car is also 
surrounded by bright, extended X-ray emission 
that appears to be assciated with the Carina Neb- 
ula: it is comparable in X-ray temperature (kT), 
size and morphology to the X-ray emission around 
HD5980 (Seward & Mitchell 1981; Fig.8). But it 
is generally assumed that this Nebula has been 
formed through the collective action of all stars of 
the Carina cluster, not only by the single action 
of n Carinae. 

We have also searched for clues of an interac- 
tion between HD 5980 and its surroundings (e.g. a 
LBV nebula). Walborn (1978) noted evidence of 
such interaction: an arc of radius ~ 30" centered 
on the star is actually visible in Ha and [O in] (but 
not in [S11]). However, this arc is rather diffuse 
(see Fig. 2), indicating a lack of shock compres- 
sion. Moreover, its velocity coincides with that 
of the quiescent region, and no line-splitting re- 
gion is seen beginning at the arc and extending 
towards the star (Danforth et al. 2002): this 
arc is certainly not the rim of an expanding shell. 
This optical arc has different size, shape, and lo- 
cation from an X-ray bright arc (Fig. 2), therefore 
we conclude that they are probably two distinct 
features unrelated to each other. Echclle spectra 
of this area also do not show any enhancement 
of the [Nil] /Ha ratio which usually indicates the 
presence of circumstellar material surrounding the 
star. The comparison between X-rays and visible 
data thus do not show any clear indication of an 
interaction between HD 5980 and its environment. 
We also note that no circumstellar nebula pro- 
duced by HD 5980 could explain the rcdshifted UV 
absorptions in the spectra of the star, or the non- 
thermal radio emission from the extended source. 

The last possibility is that the source is indeed 
a SNR, but that HD 5980 is still directly respon- 
sible for it. Several authors have proposed the 
existence of a third component in HD 5980, 'star 
C, to interpret the visual light curve (Breysacher 
& Perrier 1991) and the presence of stable pho- 
tospheric absorptions in UV (Koenigsberger et al. 
2000). To explain the existence of a SNR, an ad- 
ditional, fourth component of the system, 'star D', 
is needed. It would now be a compact object, a 
remnant of the star that exploded in SN long ago. 
If star C is much farther away than the system 



A+B+D, then the high velocity rcdshifted UV 
absorptions could be produced by star C shining 
through the receding back wall of the SNR. The 
presence of a compact object might explain the 
~6 h variations seen in the lightcurve of HD 5980 
by Sterken & Breysacher (1997). Accretion onto 
this compact object might also constitute another 
source of X-rays, thus contributing to understand 
the high X-ray emission from the system. 

6. Summary 

In this series of articles, we report the analysis 
of the Chandra data of N66, the largest star for- 
mation region of the SMC. In this first paper, we 
have focused on the most important objects of the 
field: the NGC 346 cluster and HD 5980. 

The cluster itself is relatively faint, with a total 
luminosity of L u £ ahs ~ 1.5 x 10 34 erg s" 1 in the 
0.3-10.0 keV energy range. Most of this emission 
seems correlated with the location of the brightest 
stars of the core of the cluster, but the level of X- 
ray emission probably cannot be explained solely 
by the emission from individual stars. 

In this field lies another object of interest: 
HD 5980, a remarkable star that underwent a 
LBV-type eruption in 1994. Chandra is in fact 
the first X-ray telescope to detect HD 5980 indi- 
vidually. In X-rays, the star appears very bright, 
comparable only to the brightest WR stars in the 
Galaxy, but still fainter than 77 Carinae. The com- 
parison of our results with future X-ray observa- 
tions will enable us to better understand HD 5980: 
for example, phase-locked variations will be anal- 
ysed in the perspective of the colliding wind be- 
haviour of this binary, while other variations may 
be related to the recent LBV eruption. Follow- 
up observations are thus needed to complete the 
study of this system. 

A bright, extended X-ray emission is seen to 
surround HD 5980. It is most probably due to a 
SNR whose progenitor is unknown. The spatial 
coincidence of this extended X-ray emission with 
the peculiar massive star suggest an association 
between these two objects. We also note the close 
resemblance of this X-ray emission to the Carina 
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Nebula in which the LBV rj Carinae lies. 
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Fig. 1 . — Chandra color image of the region close 
to HD 5980. Three energy bands were used to 
create this color image : red corresponds to 0.3- 
1.0 keV, green to 1.0-2.0 keV and blue to 2.0- 
10.0 keV. Before combination, the images were 
smoothed by convolution with a gaussian of a=l". 



Fig. 2.— 77STWFPC2 Ha data of the region sur- 
rounding HD 5980 (which is the star located at the 
center of the PC chip). The X-ray data overlaid as 
contours on the HST image were first binned by 
a factor of 2 and then smoothed by convolution 
with a gaussian of a=2" . 



Fig. 3. — The Chandra X-ray data of the 
NGC 346/HD 5980 region superimposed on the 
DSS image, showing the extended nature of the 
X-ray source associated with the cluster. The lo- 
cation of the brightest stars in the cluster are in- 
dicated. The X-ray data were first binned by a 
factor of 2 and then smoothed by convolution with 
a gaussian of a = 2") 

Fig. 4.— The X-ray spectrum of the NGC 346 
cluster, shown along with the best-fit absorbed 
mekal model, with N(H) = 0.42 x 10 22 cm~ 2 and 
kT = 1.01 keV (see § 3). 

Fig. 5. — a. The Chandra X-ray spectrum of 
HD 5980, shown along with the best-fit absorbed 
mekal model, with N(H) = 0.22 x 10 22 cm~ 2 and 
kT = 7.04 keV. b. The Chandra X-ray lightcurve 
of HD5980, in the 0.3-10 keV range and with 10 
bins of 10 ks each. 

Fig. 6. — The Chandra X-ray spectrum of the 
extended emission around HD 5980, shown along 
with the best-fit absorbed mekal model, with 
N{H) = 0.12 x 10 22 cm" 2 , kT = 0.66 keV and 
Z = O.UZq (see § 5). 



Fig. 8. — X-ray images of the Carina Nebula 
(24 ks ; ROSAT PSPC observation #900176) com- 
pared to the extended emission around HD 5980. 
If we assume a distance of 2.3 kpc for the Ca- 
rina Nebula and 59 kpc for the SMC, both images 
are ~55 pcx55 pc. ROSAT PSPC support rings 
and Chandra ACIS-I CCD gaps can be spotted 
on these images. Both images were smoothed by 
convolution with a gaussian (er=15" for the Ca- 
rina Nebula and a=2" after binning by a factor of 
2 for HD5980). 



Fig. 7. — Narrow-band X-ray images of the ex- 
tended emission surrounding HD5980. Each im- 
age is labeled with the principal line-emitting ion 
and the energy range used to generate the image. 
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